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Abstract — In this work, we identify a physical mechanism 
for the phenomenon of negative output differential resistance 
(NDR) in short-channel graphene FETs (GFETs) through non- 
equilibrium Green's function (NEGF) simulations. This NDR 
phenomenon is due to a transport mode bottleneck effect induced 
by the source and drain contacts, which we further elucidate 
using a simpler semianalytical ballistic model that captures the 
essential physics. This NDR effect is found to occur only when the 
gate biasing produces an n-p-n or p-n-p polarity configuration 
along the channel, for both positive and negative drain-source 
voltage sweep. In addition, we also explore the impact on the 
NDR effect of contact-induced energy broadening in the source 
and drain regions and a finite contact resistance. 

I. Introduction 

GRAPHENE has attracted considerable interest in recent 
years for applications in analog radio frequency (RF) 
electronics GQ-|4|. The reason lies in the fact that the high 
carrier mobility and Fermi velocity of graphene could allow 
device operation up to the THz range of frequencies, while 
the small on-off current ratio resulting from the zero bandgap, 
which currently prevents the use of graphene in digital elec- 
tronics, does not pose a problem in principle for analog appli- 
cations . Indeed, an integrated RF circuit made of graphene 
devices has already been demonstrated [6|. However, in more 
general analog circuits, devices with current saturation, i.e. 
small output conductance gy, are usually required. This is 
because the intrinsic voltage gain g m /g<i, where g m is the 
device transconductance, must be large. Unfortunately, current 
saturation in graphene devices is not easily obtained due to the 
lack of a bandgap. 

A quasi-saturation of the output characteristics has actually 
been reported for some experimental long-channel devices Q- 
|9l . This quasi-saturation is commonly attributed to a charge 
"pinch-off" effect due to the crossing of the quasi-Fermi level 
with the channel potential Q, ifTol — lfl~3l . Recent experiments 
have shown that not just quasi-saturation (gd — > 0) but also 
negative differential resistance (gd < 0) is possible lfl4l . lfT31l . 
Indeed, the NDR effect can also be explained with a simple 
charge "pinch-off" argument |14| within a diffusive transport 
framework, where the quasi-saturation arises as a particular 
case. Hence, understanding the origin of NDR is crucial to 
the engineering of gd in graphene. 
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In this work, we identify a different physical mechanism for 
NDR due to the role of contacts, for short-channel graphene 
transistors operating in the ballistic regime. Indeed contacts 
lfl6ll were recently found to dominate transport properties of 
graphene transistors in the ballistic regime. Although there 
have been quantum transport studies of the NDR effect in 
the ballistic regime |[T7l - |[T9l . its origin was not completely 
elucidated, in particular with regard to the effect of contacts 
and self-consistent electrostatics, which we found to play an 
important role in the operating regimes of NDR. It is the 
purpose of this worlfl to clarify the origin of NDR in ballistic 
GFETs and provide guidance to future experiments. 

The paper is organized as follows. Section [TT] reviews 
the current interpretations of the quasi-saturation in long- 
and short-channel GFETs, which prepares the ground for 
understanding the mechanism behind NDR. Then, Section [HI] 
describes our simulation models for GFETs. The results are 
shown in Section IIV1 followed by a discussion in Section [V] 
Conclusions are finally drawn in Section [VTl 

II. Review of quasi-saturation and NDR in long- 

VS. SHORT-CHANNEL GFETS 

For long-channel GFETs, the phenomenon of quasi- 
saturation and NDR can be explained using drift-diffusion 
models, which assume semiclassical diffusive transport Q, 
iflOl . ifTTl . The physical picture is sketched in FigQJ: electrons 
and holes are assumed to be described by a single quasi-Fermi 
level Ep(x), a function of the position x along the channel; 
assuming a thin gate oxide, the position in energy of the 
Dirac point in the channel region Ed follows the gate-source 
voltage Vqs ( we neglect the fact that also Ed varies with 
position), while the drain-source voltage Vds is responsible 
for the bending of Ep(x) from the Fermi level in the source 
Us to the Fermi level in the drain fip>. The drain current I is 
given by 

I = q[n{x)+p(x)]v(x)W, (1) 

where q is the electronic charge, n(x) and p(x) the sheet 
concentration of electrons and holes respectively, v(x) the 
common drift velocity along the transport direction, and W 

'This paper is an extended version of |20|. It contains new results pertaining 
to the possibility of NDR in the negative Vds regime, the effect of contact- 
induced energy broadening 1 19 1 and contact resistance, the latter being likely 
the reason why experimental short-channel devices show only weak NDR 
GS. 
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Fig. 1. Band profile in a GFET in the diffusive (a) and ballistic (b) regimes. 
The cones represent the graphene dispersion relation E(k) = Ed ± hvp\k\ 
and are used to indicate the position in energy of the Dirac point Ed- Also 
shown in (b) is the plot of M(E). In both the diffusive and ballistic cases, 
quasi-saturation occurs when the drain Fermi level crosses the channel Dirac 
point. 



the channel width. The output conductance then consists of 
two contributions: 
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While v(x) generally increases with Vbs due to the increasing 
bending of the quasi-Fermi level, the sum [n(x) + p(x)} 
instead is decreasing when Ep{x) approaches and eventually 
crosses Ed at the drain side. This charge "pinch-off is the 
effect used to explain the quasi-saturation in the diffusive 
regime. It is also possible that velocity saturation, due to 
scattering with substrate polar phonons, contribute to quasi- 
saturation [7 1, [21], since it implies dv(x) /<9Vds — > in the 
second term in the r.h.s. of Eq. [2] However, there is no clear 
consensus regarding the velocity saturation effect at present 
ll8l . In addition, the sign of gd as given by Eq. [2] can actually 
be negative if the first term prevails over the second. Indeed, 
the same arguments used for explaining quasi-saturation in 
long-channel devices can also explain NDR lfl4l . 

For ballistic GFETs instead, quasi-saturation can be under- 
stood using the Landauer formalism lfl2l (Fig. QJ). Assuming 
for simplicity the zero-temperature approximation, the energy 
window for transport is the one between (ip and /is and the 
current is given by 



2q f^ s 
J = -M M(E)dE, 



(3) 



where h is Planck's constant, M(E) = 2W \E - E d \ /{irhv F ) 
the number of propagating modes in the graphene channel at 
energy E, h — h/(2n), and vp the graphene Fermi velocity. 
In this case, the output conductance is 

dl .„ , dM(E) 
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M{E) is shifted up or down in energy by varying Ed, which 
in turn is determined by the gate electrostatics. Assuming that 




Fig. 2. Semianalytical model: Dirac point profile (top) and population of 
states in fc-space (bottom) corresponding to the indicated energy E. States 
populated with the Fermi function of the source/drain are indicated with a 
thick red/blue (light-gray /dark-gray) line. 



Ed does not depend on Vds> the second term in the r.h.s. of 
Eq. 2] is zero. Consequently, it can be seen that the current 
tends to saturate when [ip, approaches the Dirac point in the 
channel because M(Ed) = (here we neglect the graphene 
minimum conductivity |22|, [23|). The sign of m Eq. |4] 
can be negative only if dM(E)/dVos is negative, which is 
not obvious from this simple picture. In the following, we will 
generalize this simple model to include the effects of contacts 
and self-consistent electrostatics, and show that dM(E) / <9Vds 
can indeed be negative within the energy window for transport. 

III. Models 

In this section, we describe in detail the semianalytical 
model and the numerical model which are used to investigate 
the device behavior. 

The semianalytical model assumes a simple ideal square 
potential barrier, as shown in Fig. ID-top, where Ed y s/D/G 
is the energy of the Dirac point in the source/drain/channel 
region. This approximation neglects the tunneling probability 
effect due to the finite transition length of the step potential 
[24 1, [123, which is accounted for in the numerical treatment 
later. Metal induced doping |26| is introduced through a fixed 
difference AE con between the Fermi level and the Dirac point 
in the metal-doped source and drain regions (Fig. [2]- top). 

The transport model is composed of the following equations 
for the electron and hole concentration in the channel region 
and for the drain current: 



// = / dE [D L (E)f s (E) + D R (E)f D (E)} , 

rE d , G 

V = ' 



(5) 



dE [D L (E) (1 - fs(E)) + Dji(E) (1 - f D (E))] , 

-oo 

I=jj dEM(E)[f s (E)-f D (E)}, 



(6) 
(7) 



where Dl/r(E) is the density of states (DOS) in the chan- 
nel at energy E relative to injection from source/drain and 
Is/d{E) the contact Fermi distribution with Fermi level hs/d- 
In turn, the model for D L (E), D R (E), and M(E) is given 
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where d is the Heaviside step function and 



-Os)] 



M 



S/D/G 



(E) = 



2W 



S/D/G ■ 



(8) 

(9) 
(10) 

(ID 



The numerical value for vp is set using the equation vp = 
(3/2)acc|^|/^ in which acc = 1-42 A is the carbon-carbon 
distance in graphene and t = —2.7 eV the tight-binding 
parameter describing hopping between nearest neighbor p z 
orbitals. The quantities ks/D/G an d @S/d are defined as (their 
physical meaning is illustrated in Fig. |2]-bottom) 



E-E, 
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The Dirac point in the channel E dG is self-consistently 
computed with n and p through a plane-capacitor model which 
accounts for electrostatics: 

q(n -p) = C ox {-us jq + V GS + E d . G /q) , (14) 

where C ox is the gate oxide capacitance and a zero workfunc- 
tion difference is assumed between gate and graphene. 

The model in Eqs. [514X31 corresponds to the solution of the 
ballistic Boltzmann equation in the channel region assuming 
energy and transverse momentum conservation at the two 
junctions. As an example, Fig.|2]-bottom shows the distribution 
function in fc-space corresponding to the potential in Fig.[5]-top 
and at the indicated energy E. The red/blue (light-gray/dark- 
gray) color represents f s / D (E). The plot can be understood 
by assuming that the transmission probability across each 
junction, for an incident electron with transverse momentum 
k y , is either 1, if states with the same k y are available on 
the other side of the junction, or otherwise. In the figure, 
electrons from the source (red or light gray) are perfectly trans- 
mitted through both junctions, thus populating only rightward 
propagating states in the channel (note that the group velocity 
is opposite to k for states in the valence band). The ones from 
the drain (blue or dark gray) enter the channel with probability 
one; at the source-channel junction, they are either perfectly 
transmitted if \k y \ < k$ or totally reflected if \k y \ > k$, thus 
populating both left- and rightward propagating states. 

We highlight the fact that the current contribution at a given 
energy is determined by the region where the Fermi surface 
has the smallest radius (Eq. [TOl . which means a transport 
bottleneck effect due to the series of graphene junctions. This 
model for M{E) was first discussed in |24| and was also used 
to describe NDR in single p-n junction devices l27l . It is worth 
noting that, if ks,kp> > k G (i.e., if the number of modes in 
the contacts is larger than in the channel), the same model for 
charge and current as in lfl2l is recovered. 



To benchmark the semianalytical model, we use an atomistic 
full-quantum code for graphene nanoribbons |28|, based on 
the self-consistent solution of the tight-binding (TB) NEGF 
and 3D Poisson equations and optionally including graphene 
acoustic phonon (AP) and optical phonon (OP) scattering. The 
source and drain regions are treated as in the semianalytical 
model with a fixed AE'con and semi-infinite extensions (as 
in Ull). In the following, we will show results from both 
ballistic simulations and simulations with phonon scattering; 
in the latter case, we use the parameters -D^.ap = 0.03 
eV 2 , D U:0 p = 0.027 eV 2 , and huj p = 160 meV, whose 
definitions can be found in E9\ . 

In the semianalytical model presented above, we have 
assumed that the source and drain regions are described by 
the same conical electronic dispersion relation as the channel 
(Eq.fT2l. In reality, the graphene DOS in the contacted regions 
is broadened due to the coupling with the metal contacts, so 
that a finite DOS (and thus a finite current injection) is induced 
at the Dirac point. In the following, we study separately 
the effect of contact-induced energy broadening. Regarding 
the NEGF code, we include the broadening as a constant 
imaginary diagonal self-energy — iA for the source and drain 
regions ||T9"1 . We have verified that the resulting DOS in the 
source/drain region, Dgi D (E), can be well reproduced by the 
formula 



v (E-E diS/D ) + A 2 

D S /d(E) = 2^ — (15) 

7T (nvF) 

where A is a fitting parameter. Assuming the same relation 
between Mg/p, and Ds/d as in tne case without broadening, 

M S/D (E) = D s/D (E)hv F W, (16) 

we get an effective dispersion relation 



k s/D = —M S/D (E) 



(E-E, 



d,S/D) 



A 2 



flVl 
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which we use in place of Eq.[T2]for k s /p> to capture the effect 
of energy broadening within the semianalytical model. 

IV. Results 

We consider n-type doped source and drain regions (unless 
stated otherwise, the values AE'con = 0.4 eV and A = A = 
are assumed) and an ideal gate dieletric with equivalent oxide 
thickness (EOT) of 0.5 nm. All the results presented in the 
following are at room temperature. With Vqs = and at 
equilibrium, the channel Dirac point is aligned with the source 
and drain Fermi levels and the channel is intrinsic; by applying 
a positive (negative) Vgs, the bands in the channel are shifted 
down (up) thus creating an n-n-n (n-p-n) double junction. 
We explore in the following both the n-n-n and n-p-n bias 
conditions. 

A. Quasi-saturation in n-n-n structure with Vds > 

In Fig. [3] we plot the output characteristics for Vqs > 
(n-type channel) from the semianalytical model and the one 
obtained by setting ks,kp> > k G and thus 9$ = Op, = tt/2 
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Fig. 3. Normalized output characteristics for Vgs > (n-n-n configuration) 
from the semianalytical model (solid line) and the one obtained by assuming 
^Si^-D > (dashed line). 
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Fig. 4. Normalized output characteristics for Vgs < (n-p-n configuration) 
from the semianalytical model (solid line) and the one obtained by assuming 
^Si^-D > kc (dashed line). 
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Fig. 5. Dirac point profile (left), number of propagating modes versus energy 
(center), and current spectrum (right) from the semianalytical model, for two 
different Vbs at Vgs = -0.4 V. 




Fig. 6. Normalized output characteristics for Vgs < from the semianalyt- 
ical model and ballistic NEGF simulations of nanoribbons with two different 
widths W . 



in Eqs. [8UTU1 (cfr. [12]). It can be seen that the two models 
give similar results at large Vgs> both predicting the quasi- 
saturation behavior discussed above. This means that the mode 
reflection at the two junctions has no significant effect in this 
bias condition. However, at lower Vgs> the two models depart 
significantly for large Vbs- Indeed, at small Vgs an d large 
Vbs> the channel doping is actually turned p-type by the drain 
contact and the transport regime is similar to the n-p-n case 
discussed below. 

B. NDR in n-p-n structure with Vbs > 

The output characteristics for Vgs < (p-type channel) are 
shown in Fig. [4] While the model neglecting the finite number 
of modes in the source and drain predicts a monotonically 
increasing current, the one proposed here clearly gives NDR. 

The origin of the NDR effect is explained by looking at 
Fig- El which compares the band profile, number of modes, and 
current spectrum (integrand in Eq. [7]l obtained at two different 
Vbs biases along the Vgs = —0.4 V curve. NDR is the combi- 
nation of two effects. First, for energies close to E^ g, M{E) 
is limited by the number of modes in the source: in particular, 
for E — Ed,s we nave M(E) = 0, and this leads to the quasi- 
saturation behavior of the current for jio approaching Ed,s, 



as already pointed out in [19] and mathematically represented 
(in the zero-temperature approximation) by the first term in 
the r.h.s. of Eq. 0] Secondly, by lowering /zd, the flux of 
electrons injected from the drain into the channel is reduced, 
while the flux from the source is kept fixed, causing a hole 
pile-up in the channel: due to the electrostatic feedback, Ed,G 
is lowered and M(E) decreases for a portion of the energy 
range within /ig and hd, leading to a decrease of the current 
rather than a saturation, as expressed (in the zero-temperature 
approximation) by the second term in the r.h.s of Eq. |4] The 
fact that the decrease in M{E) is not fully compensated by 
the larger (/xg — fijj) is confirmed by the plot in Fig. |5}right, 
where the area under the red curve (current at higher Vbs) is 
slightly smaller than the area under the black curve (current 
at lower Vbs)- 

We highlight the necessity of two Dirac points, one in the 
channel and the other one in either the source or the drain, 
for NDR to be possible instead of saturation. Also, we note 
that a similar explanation for NDR was given in ifTHl . even if 
a simple, not self-consistent model for barrier lowering was 
used (shift of Ed,G by — Vbs/2 with respect to the value at 
Vbs = 0). 

The qualitative shape of the output characteristics is con- 
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Fig. 7. Dirac point profile (left) and current spectrum (right) from the 
semianalytical model and from ballistic NEGF simulations for the same two 
biases as in Fig. \5\ 
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Fig. 8. Normalized output characteristics for Vqs ^ from NEGF 
simulations with ballistic transport, with AP scattering, and with both AP 
and OP scattering. 
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Fig. 9. Normalized output characteristics for Vgs < an d for both 
positive and negative Vbs from the semianalytical model and ballistic NEGF 
simulations. 
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Fig. 10. Dirac point profile (left), number of propagating modes versus 
energy (center), and current spectrum (right) from the semianalytical model, 
for two different negative Vds at ^GS = —0.4 V. 



firmed by ballistic NEGF simulations (Fig. [6) considering a 
50-nm-channel length. Since the NEGF code we use is for 
graphene nanoribbons, we have checked the convergence of 
the results with respect to the nanoribbon width 147; the two 
sets of NEGF curves for W — 20 nm and W = 30 nm in Fig. [6] 
overlap, indicating that the nanoribbons are wide enough to 
behave as 2D graphene (in real devices, larger widths would 
be needed to reach the 2D behavior because of the effect of 
edge roughness). 

In Fig. [7] we plot the band profiles and current spectra 
obtained from the two models. In the ballistic NEGF formal- 
ism, the number of propagating modes M(E) in the current 
expression is replaced by the transmission function T(E) = 
^2j—i Tj(E), where Tj(E) is the transmission probability 
from source to drain of mode j |30|. From the figure, it can 
be seen that the assumption of square potential barrier is well 
justified and that the barrier lowering is similar; the lower 
current spectrum in the NEGF case can be explained with 
wavefunction mismatch at the junctions, causing Tj(E) < 1 
for propagating modes even in the case of an abrupt potential 
step El, l3fl . 

Finally, we have investigated the effect of scattering due to 
graphene longitudinal acoustic and optical phonons and found 



that, at this channel length, it is too weak to affect the current 
(Fig. |8j. Nevertheless, it is worth noting that the effect of OP 
scattering is to increase slightly the current in the NDR region 
with respect to the case with only AP scattering: this can be 
explained by the fact that OP scattering creates additional 
paths for the current flow due to energy relaxation, thus 
reducing the transport bottleneck effect induced by M(E). 

C. NDR in n-p-n structure with Vds < 

We have also studied the device behavior when the polarity 
of the drain voltage is reversed. Fig. [9] shows the complete 
output characteristics for Vgs < including both positive and 
negative Vds values (the curves for Vds > are the same 
as those shown in Figs. [4] or [6j. Another NDR effect, not 
previously reported, is observed at Vds < in both the results 
of the semianalytical model and NEGF simulations. 

To help understand the origin of the phenomenon, we plot 
in Fig. [10] the band profile and spectra for two negative Vds 
values along the Vgs = —0.4 V curve. The mode bottleneck 
induced by the Dirac point at E = E^g is responsible for 
the current saturation as fijj approaches E^ q, since it gives 
M (hd) in Eq. [4] At the same time, when fijj is raised, 
the mode bottleneck at E = Ed,D causes a decrease of M (E) 
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Fig. 11. Dependence on Vgs and Ai? C on of the voltage swing (left) and 
peak-to-valley current ratio (right) of NDR at Vbs < 0. The results are from 
the semianalytical model. The inset shows the definition of V ma x, V m i n , 
-/max, and Xmin* 
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Fig. 12. Normalized output characteristics for Vgs < from the semiana- 
lytical model and ballistic NEGF simulations including energy broadening in 
the source and drain. 

within the energy window between fi$ and //£>, due to the rigid 
shift of Ed.D with /j,rj- This corresponds to dM (E)/dVas < 
in Eq. [4] so that NDR rather than saturation occurs. 

We conclude that the mechanisms for NDR at the two drain 
voltage polarities are very similar and essentially related to the 
combined effect of two Dirac points, one in the channel and the 
other in either the source or the drain, in limiting the current. 
However, NDR at Vbs < does not involve an electrostatic 
feedback in the channel and it is thus expected to be a more 
robust phenomenon than its counterpart at Vbs > 0. 

While the optimization of NDR is beyond the scope of 
this work, we just show in Fig. QT| how NDR at Vbs < 
can be modulated by varying the electrostatic doping of the 
contacts, i.e. AE con . From the figure, it can be seen that 
both the voltage swing |y max — V m - m \ and the peak-to-valley 
current ratio imax/Jmin (the symbols are defined in the inset 
of Fig. fTTI) can be somewhat enhanced by increasing the 
contact doping, which also enlarges the Vgs window where 
the phenomenon occurs. 

D. Effect of energy broadening in the source and drain regions 

To study the effect of energy broadening due to metal- 
graphene coupling in the source and drain regions, we set 



Fig. 13. Normalized output characteristics from the semianalytical model 
for different values of the contact resistance Rc, at Vgs = —0.4 V (left) 
and Vqs = — 0.6 V (right). Rq is split equally between the source and drain 
terminals as shown in the inset. 

A = 50 meV in the NEGF code and A = 162 meV in the 
semianalytical model. The latter has been fitted to provide 
the same value of DOS at the Dirac point, Ds/d(Eci,s/d) ~ 
3 x 10 13 cm~ 2 eV _1 , as in the NEGF case. The resulting output 
characteristics are plotted in Fig. Q~2] The two methods give 
qualitatively similar results. At Vbs > 0, NDR disappears as 
already observed in |[T9l . However, we find here that NDR, 
albeit weak, is still possible at Vbs < 0. The experimental 
verification of the effect should thus be easier in the Vds < 
case. 

E. Effect of contact resistance 

We finally consider the effect of a series contact resistance 
Rc equally split between the source and drain terminals (inset 
of Fig. [T3]l. The semianalytical model is modified accordingly, 
by replacing the Fermi level [is ID entering the equations for 
fs/D and E,i s/d with the Fermi level ^5/ of the intrinsic 
source/drain. The latter is calculated self-consistently with the 
current / through the resistor equation 

Rc 

fJ-S'/D' = fJ-S/D T I — 1 ' ( 18 ) 

where the upper/lower sign holds for source/drain. The / vs. 
lbs characteristics for different values of Rc and two values 
of Vqs are shown in Fig. [13] It can be seen that NDR tends 
to disappear with increasing Rc due to the increasing voltage 
drop across the two resistors. Already at Rc = 300 • /Ltm, 
which is a typical experimental value [14], the output curve 
resembles the linear characteristic of a resistor, for both 
values of the gate voltage considered, indicating that contact 
resistance is a major problem for the operation of short- 
channel devices. 

V. Discussion 

The analysis above indicates that, given a specific type of 
doping of the contacted graphene regions, the gate voltage 
needs to be biased such that an n-p-n or p-n-p double junction 
is formed, in order for NDR in ballistic GFETs to be possible. 
In addition, NDR is expected to be more likely for drain 
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voltages with the same polarity as Vbs> due to the higher 
robustness of NDR against energy broadening in this case. It 
is interesting to compare these findings with the actual bias 
conditions employed in experiments to observe NDR. 

We note that the long-channel devices in |14|, 031 were 
biased in an n-n-n or p-p-p configuration thus ruling out our 
interpretation of the phenomenon. Indeed, for channel lengths 
of the order of 1/im, scattering is expected to cause significant 
momentum and energy relaxation and charge inhomogeneity 
along the channel, so that the drift-diffusion interpretation 
discussed in Section [TT] seems more appropriate. 

On the other hand, in the Supporting Information of fl4l . 
results for GFETs with channel lengths of about 100 nm have 
been reported too. Some of these devices show NDR for gate 
voltages around V: this biasing could be the counterpart case 
of the one mentioned at the end of Section llV-AI (incompletely 
formed p-type channel, p-type contacts, and Vds < 0), 
so that our explanation of NDR could apply in this case. 
However, further experimental evidence is needed to prove 
this hypothesis (for example by exploring all the Vgs an d 
Vds polarities). 

VI. Conclusions 

Through a semianalytical model and detailed quantum sim- 
ulations, we have clarified the nature and bias conditions for 
NDR in short-channel GFETs. The origin of the phenomenon 
is attributed to the transport-mode bottleneck induced by the 
graphene Dirac point. The combined effect of two Dirac 
points, one in the channel and the other in either the source 
or the drain, is necessary for NDR to occur instead of quasi- 
saturation. This is verified in the n-p-n or p-n-p configuration, 
for both polarities of Vds- In the presence of energy broad- 
ening due to the metal-graphene coupling in the source and 
drain regions, NDR disappears at one Vds polarity, but a weak 
NDR, or at least saturation, is still attainable at the other one. 
It is also found that contact resistance at typical experimental 
values suppresses NDR, representing a major obstacle for the 
verification of the phenomenon in experiments. 

The NDR mechanism could offer new possibilities for 
the optimization of the saturation behavior of the output 
characteristics of analog GFETs. The semianalytical model 
presented here, providing a good physical insight of NDR, 
could be a useful simulation tool for such an optimization 
study. 
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